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pN . Abstract 

^ , Generalized Drinfeld-Sokolov (DS) hierarchies are constructed through local reduc- 

1^ . tions of Hamiltonian flows generated by monodromy invariants on the dual of a loop 

algebra. Following earlier work of De Groot et al, reductions based upon graded regular 
elements of arbitrary Heisenberg subalgebras are considered. We show that, in the case of 
the nontwisted loop algebra i{gln), graded regular elements exist only in those Heisenberg 
subalgebras which correspond either to the partitions of n into the sum of equal numbers 
77, = pr or to equal numbers plus one n = pr + 1. We prove that the reduction belonging 
to the grade 1 regular elements in the case n = pr yields the p x p matrix version of the 
Gelfand-Dickey r-KdV hierarchy, generalizing the scalar case p = 1 considered by DS. The 
methods of DS are utilized throughout the analysis, but formulating the reduction entirely 
within the Hamiltonian framework provided by the classical r-matrix approach leads to 
some simpliflcations even for p = 1. 
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0. Introduction 

The generalized KdV type liierarctiies of Drinfeld and Sokolov (DS) are among ttie 
most important examples in the field of integrable evolution equations [1] . They also play 
an important role in current studies of two-dimensional gravity [2] and in conformal field 
theory [3] . The "second Gelfand-Dickey" Poisson bracket of these bihamiltonian systems 
is a reduction of the affine current algebra Lie-Poisson bracket, and it gives an extension 
of the Virasoro algebra by conformal tensors. Such extended conformal algebras are called 
W-algebras and have received a lot of attention recently [4-6] . 

The motivation for the present work was to gain, from a purely Hamiltonian viewpoint, 
a better understanding of the reduction procedure used in [1] and the generalizations 
proposed in a recent series of papers [7-9] aimed at the construction of new integrable 
hierarchies and W-algebras. 

In the DS construction [1] one starts by considering a first order matrix differential 
operator of the form 

C = d + ^, with fx = q + A, (0.1) 

where A is a matrix representing a grade 1 regular element of the principal Heisenberg 
subalgebra of a loop algebra, and g is a smooth mapping from 5"^ into an appropriate 
subspace of the loop algebra represented by "lower triangular matrices" . The crucial step 
of the construction is to transform n into the Heisenberg subalgebra by a conjugation of 
C That this can be achieved by an algebraic, recursive procedure is due to the fact that 
the grades in q are lower than the grade of A and A is regular. The fact that A is regular 
also implies that the stabilizer of the transformed operator is given by the Heisenberg 
subalgebra. The compatible zero curvature equations are obtained from the positively 
graded generators of this Heisenberg subalgebra by transforming them into the stabilizer 
of C and applying a splitting procedure. The system exhibits a gauge invariance under a 
nilpotent, "strictly lower triangular" , gauge group. 

The authors of [7] realized (see also [10], [11]) that the DS construction can be applied 
in more general circumstances. They proposed to derive new hierarchies by replacing A in 
the above by any positive, graded regular element of any graded Heisenberg subalgebra, 
and correspondingly modifying the DS definition of the variable q and the gauge group. 
However, apart from some very simple examples, they did not investigate which systems 
can be obtained on the basis of this rather general proposal. We shall see here that the 
number of new hierarchies arising from this approach, which in [7] were termed type I 
hierarchies, is in fact rather limited, since graded regular elements do not exist in most 
Heisenberg subalgebras. For simplicity, we shall investigate here the case of the nontwisted 



loop algebra based on the general linear Lie algebra gin, in which case the inequivalent 
Heisenberg subalgebras are classified by the partitions on n [12]. By using the explicit 
description of the inequivalent Heisenberg subalgebras given in [13], we shall prove that 
graded regular elements exist only for the special partitions 

n = pr and n=pr + l. (0-2) 

After explaining this observation, we shall give a detailed analysis of the first series 
of "nice cases" under (0.2). The grade 1 regular elements of the corresponding Heisenberg 
subalgebra are the n x n matrices of the form Ar,p = A^ CS -D, where D is a p x p diagonal 
matrix such that D^ has distinct, non-zero eigenvalues, and A^. is the usual rxr "Drinfeld- 
Sokolov matrix" containing I's above the diagonal and the spectral parameter A in the 
lower-left corner. The "constrained manifold" of the generalized DS reduction will be 
taken to be the space of £'s of the form (0.1), where now A = Ar^p and g is a mapping 
from S^ into the block lower triangular subalgebra of gin, with p x p blocks. As in the 
p = 1 case considered by DS, the "reduced space" will be obtained by factorizing this 
constrained manifold by the group of nonabelian gauge transformations, 

C ^ gCg-^ , (0.3) 

where g is now block lower triangular, having p x p unit matrices in the diagonal blocks. 
We shall place the construction in a Hamiltonian setting from the very beginning, from 
which it will be clear that the reduced space is a bihamiltonian manifold that carries the 
commuting hierarchy of Hamiltonians provided by the monodromy invariants of C It will 
also be clear that the locality of the reduced system is guaranteed by construction. 

In order to describe the reduced system (i.e., the generalized KdV hierarchy) in terms 
of gauge invariant variables, we shall prove the following facts, which extend the p = 1 
results of [1]. First, the reduced space is the space of "matrix Lax operators" of the form 

L = (-D)-'^a'^ + uid"^-^ + ... + Ur-ld + Ur , (0.4) 

where the Ui are now smooth mappings from S^ into the space oi px p matrices. Second, 
the reduced bihamiltonian structure is that given by the two compatible (matrix) Gelfand- 
Dickey Poisson brackets; the second Poisson bracket algebra qualifies as a classical W- 
algebra. Third, the hierarchy of commuting flows is generated by the Hamiltonians given 
by integrating the componentwise residues of the fractional (including integral) powers 
of the pseudo-differential operators obtained by diagonalizing the matrix Lax operators. 



In short, the DS reduction extends in this case to yield the p x p matrix version of the 
well-known (e.g. [14]) Gelfand-Dickey r-KdV hierarchy. 

The main additional step required in computing the Hamiltonians in the px p matrix 
case consists in the diagonalization of L. Those Hamiltonians which are obtained from the 
integral powers of the diagonalized Lax operators generate bihamiltonian ladders whose 
first elements are the Hamiltonians given by integrating the diagonal components of ui, 
which are Casimirs of the first Poisson structure. The number of these bihamiltonian 
ladders (which are missing in the scalar case) is p — 1 since the integral of tr{D'^ui) is 
a Casimir of both Poisson structures. The other Hamiltonians can also be described as 
integrals of trace-residues of independent fractional powers of L, without diagonalization. 

The KdV type hierarchies based on matrix Lax operators of the type (0.4) have 
been investigated before in refs. [15-17] and more recently in [18]. In refs. [15-17] the 
matrix {—D)~^ (i-e-, the coefficient of the leading term of L) was required to be regular 
because this implies the existence of the maximal number of independent rth roots of L 
and corresponding commuting fiows. It is interesting to see this condition re-emerge here 
from requiring that the matrix A used in the reduction procedure be a regular element 
of the Heisenberg subalgebra. In these papers, the additional assumption was made that 
the diagonal part of ui vanishes. Setting [wi]diag = is consistent with the equations of 
the hierarchy resulting from the DS reduction and in fact corresponds to an additional 
Hamiltonian symmetry reduction. (See Remarks 2.5 - 2.7.) 

The recent preprint [18] deals with the hierarchy defined by the fractional powers 
of "covariant Lax operators", which are equivalent to operators of the form (0.4). More 
precisely, the case of the Lie algebra sin was considered, which in our approach corresponds 
to imposing the constraints tr{D'^ui) = 0. However, instead of taking a regular matrix for 
{—D)~^, the unit matrix was used. As a result, the hierarchy obtained is much smaller 
than the one following from the DS reduction using a regular element. 

To emphasize the link between the approach used in [1], [7] and that based upon the 
Adler-Kostant-Symes (AKS) construction as presented in the present work, we give the 
following elementary lemma. 

Lemma 0.1. Let A be a Lie algebra. Let hq G A* be given and X G cent (stab (/Uq)) C A 
be an element in the center of its stabilizer. Then there exists an element ip G T{A*) 
of the ring of a-d* -invariant functions on A* such that Vt/?!;^,, = X , where, as usual, an 
identification has been made between the cotangent space to A* at fxo o'^^ (A*)* ~ A. 
Conversely, if ^ & X{A*), then Vv?|^o '■= X E cent(stab(/Uo)). 

The point of this lemma is that, when applied to the loop algebra A = i{Q) of a 



Lie algebra Q, with the sphtting A = A+ + A- into the sum of subalgebras consisting of 
positive and negative powers in the loop parameter A, it implies that, starting from any X G 
cent(stab(/Uo)), with hq G (A-)* ~ {A*) + , the flow induced in (A-)* by exponentiation 
and factorization in the group is given by integration of a Lax type equation 



'dt 



±(ad*(Vy:U)±)(^) (0.5) 



with (f G I{A*). Taking the Lie algebra Q itself as a centrally extended loop algebra in 
the space variable x, eq. (0.5) becomes a zero-curvature (Zakharov-Shabat) equation. The 
commutativity of such flows is part of the AKS theorem. 

The above lemma underlies the equivalence between the DS approach [1], [7], which 
is based essentially upon cent (stab (/u)), and the AKS approach, based on X{A*). This 
equivalence is certainly known to specialists, but in this paper it will be taken as the 
starting point and all results will be derived from the AKS Hamiltonian point of view. 

This paper is organized as follows. In Sec. 1 we collect results that are relevant 
for understanding the DS type construction of compatible zero curvature equations in 
a Hamiltonian setting. In particular, in Sec. 1.1 we recall the relevant aspects of the 
AKS construction of commuting flows. In Sec. 2.2 we discuss a sufficient condition that 
can be used to obtain local ad*-invariant Hamiltonians from the asymptotic expansion of 
the monodromy matrix of an appropriate first order matrix differential operator. These 
two sections naturally lead us to look for the graded regular elements of the Heisenberg 
subalgebras of a loop algebra as the starting point for obtaining generalized DS hierarchies. 
The solution to this problem is given in the case of the nontwisted loop algebra (-{gin) by 
Theorem 1.7 in Sec. 1.3. Sec. 2 is devoted to the generalized DS reduction yielding the 
matrix Gelfand-Dickey hierarchy. The description of the reduced space is established 
in Sec. 2.1, the Hamiltonian structures are described in Sec. 2.2 and the Hamiltonians 
themselves are given in Sec. 2.3. The Poisson brackets and the first few Hamiltonians are 
computed explicitly for an example in Sec. 2.4. The main results of Sec. 2 are Theorem 
2.4 which identifies the reduced Poisson structures as the first and second Gelfand-Dickey 
Poisson structures, and Corollary 2.11 of Theorem 2.10 which gives the generating set of 
commuting Hamiltonians. The paper concludes with remarks relating the matrix Gelfand- 
Dickey hierarchies to nonabelian "conformal" and "affine" Toda systems, comments on 
W-algebras and on the case n = pr -\-l, and further remarks concerning the literature and 
some open problems. 



1. The AKS construction and local reductions 

In the first two sections we review some well-known results about the AKS approach 
in loop algebras. We shall naturally be led to considering the problem of finding all the 
graded regular elements of the Heisenberg subalgebras of gin ® C[A, A~^], which are given 
in Section 2.3. 

1.1. The AKS construction 

Here we summarize those points of the AKS (or r-matrix) approach which we shall 
need. Readers unfamiliar with the construction could consult, for example, refs. [19-21] 
for further details. 

Let ^ be a Lie algebra with Lie bracket [ , ]. Suppose that R is a classical r-matrix; 
that is, R G End^ and the bracket [A, Y]r : Ax A ^' A given by 

[X,Y]n = l[RX,Y] + ^[X,RY] (1.1) 

is also a Lie bracket. For a function (p in C°"{A*) we define its differential VqI/? G .4. at a 
point q; G vA* by 

j^^{a + t(3)\t=o = {(3,V^^) W(3eA\ (1.2) 

where ( , ) is the dual pairing. Let 

I{A*) = {^e C^{A*)\0= ((ad*A)(a), Vay:) = (a, [V^v?, A]) VA G ^ }, (1.3) 

be the set of ad*-invariant functions on A* . Note that in many cases we can think in terms 
of the Ad*-action on A* of a Lie group G corresponding to A, in which case 

I{A*) = {ipeC'^{A*)\ipiAdla) = ip{a) V(7 G G }. (1.4) 

Here ad* means the action dual to the original Lie bracket [ , ] on ^; to refer to the action 
of A on A* dual to that given by the bracket [ , ]r, we write ad|j. 

A surprisingly large number of integrable systems arise as consequences of the following 
result [20] , which is the r-matrix version of the AKS theorem, and whose proof is a direct 
application of (1.1) and (1.3). 

Proposition 1.1. The elements of I (A*) are an involutive family in C°°{A*) with respect 
to the R Lie-Poisson bracket on A* . That is, if (^^i/j E I{A*) then 

(V?, ipJRia) = (a, [Va</?, VaV]i?) = . (1.5) 
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The dynamical equation generated by the Hamiltonian </? G I {A*) through the R Lie-Poisson 
bracket has the generalized Lax form 

d=(ad*^i?V«v')(«)- (1-6) 

Consider now the special case that will be of interest in what follows (see [21]). Let 
^ be a Lie algebra and set A = i{g) = G ® C[A, A'^] = {Y.l=r^i^'\^^ ^ ^}- Let 
A+ = g®C[X], A- =^(g)A"^C[A"^], so that A = A+ + A-. We let P± be the projection 
operators defined by this splitting and set R = P+ — P_. For any r] G C[A, A~^] we define 
?7 : ^ -^ ^ by {fiX){X) = 7]{X)X{X). Then 

Rr,:=Rof] (1.7) 

defines a classical r-matrix for any i], and the corresponding Lie-Poisson brackets for dif- 
ferent r] are compatible. 

We can identify i{g)* with i{Q*) = G* ® C[X, X~^] by means of the pairing ( , ) given 
by 

(a,X):=(a(X))|_i, i.e. (5^a,A\ J^X, V) := J^ a,{X,) . (1.8) 

i j i+j = -l 

Remark 1.1. An equivalent procedure would be to keep the r-matrix fixed and define a 
dual pairing for each r] G C[A, A~^] by inserting a factor r] into the definition (1.8). In 
view of this, we see that in this case there is no essential difference between the r-matrix 
formulation and the original AKS splitting procedure. In what follows, the collection of 
commuting Hamiltonian systems determined by the elements of I{l{Q*)), together with 
the compatible Poisson brackets { , }i? will be referred to as the AKS system. 

The ad*-action of 1{Q) on 1{Q*) is given by pointwise evaluation in A. This implies 
that one can in general obtain elements of I{l{Q*)) in the following way. Let </? G I{Q*) be 
an invariant polynomial; choose any element q G C[A, A~^]. Then ipg G C°°(£(^*)), given 
by 

^,{a) := ig{XMa{X))) |_i , « = ^ a.V G £(^*), (1.9) 

i 

is in I{i{Q*)). If I{Q*) is generated by polynomials (/? then I{£{Q*)) is generated by the 
corresponding functions (fg. The Hamiltonian vector fields determined by these invariant 
functions on £{Q*) obey the following relations, 

{/,V'.xK = {/,V'.Kx' V/gC-(£(^*)), Vr?,^,xGC[A,A-i]. (1.10) 
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Notice [21] that the space 

n 

M-m,n-={ ^ UiX^\ui e g* } C i{Q*) where m > 0, n > -1, (l-H) 

is a Poisson subspace for any of the r-matrices Rk := R o fjk, where ?7a;(A) := A'^, if 
—m < k < n + 1. Furthermore, if k < n then Un is a Casimir element; i.e., constant 
under any Hamiltonian flow of the Rk Lie-Poisson bracket. Hence the afflne subspace 
J\A-m,n C A4-m,n haviug Un fixed has n + m + 1 compatible Poisson brackets, and the 
restriction of the elements of I{l{Q*)) provides a commuting family of Hamiltonians with 
respect to any of them. Note that the same computation which shows Ai-jn,n is a Poisson 
subspace of £{Q*) justifies restricting ourselves to £{Q*), which is strictly speaking only a 
subspace of i{Q)*. In both cases one simply has to check that an arbitrary Hamiltonian 
flow determined by 

d = (ad^^V«if)(a) (1.12) 

does not leave the space. 

From now on we take Q to be gin , the central extension of the algebra of smooth 
loops in gin, i.e., Q = { {X, a) \ X : S^ ^' gin , a G C } with* 

[{X,a),{Y,b)]= (xY -YX , f dxtrX'{x)Y{x)j . (1.13) 

As usual, we represent (a dense subspace of) the dual space Q* as the space of flrst order 
matrix differential operators 

C={ed + ^)< — >{^,e)eg*, (1.14) 

with n G gin ~ gin , e G C and dual pairing 

((/U, e) , (X, a)) = ea + / dxtr ij,{x)X{x) . (1-15) 

Jo 

Upon introducing the spectral parameter A, the elements of £{Q) and £{Q*) are given by 
pairs (X^a) and (/U,e), where X = ^XjA-', /U = J2 f^i^^ ^^^ ^^^ mappings from 5"^ into 
i{gln), and a = ^ ajX^ , ^ = Yl ^i^^ ^^^ elements of C[A, A~^]. The ad*-action on £{Q*) is 
given by 

ad*(X, a) : {^, e) ^ {Xfi - fxX - eX', 0) . (1.16) 



* The periodic space variable parametrizing 5"^ is usually denoted by x G [0, 27r], and 
tilde signifies here "loops in x" . 



The Ad*-action is given by "nonabelian gauge transformations" ; that is, by 

C^ed + ^i^ gCg-^ = ed + {g^ig-^ - eg'g-') , (1.17) 

where g is an element of the loop group i{Gln) associated to i{gl^). 

The ring of ad*-invariant functions I{£{Q*)) can be specified as follows. Consider the 
linear problem 

(ea^ + ;u(x,A))$(x,A) = 0, (1.18) 

where $(x, A) G Gin- Then the eigenvalues of the monodromy matrix 

T(A)=$(27r,A)($(0,A))-\ (1.19) 

viewed as functions on i{G*), generate the ring I{i{Q*)). 

The number of invariants is infinite due to the parametric dependence of T on A. In 
general these are not local functionals of //; i.e., they cannot be expressed as functions of a 
finite number of integrals of local densities in the components Hi{x) of /u(x. A) = J2i IJ'i{-c)X^ 
and their derivatives. However, under certain conditions (see Sec. 1.2), it is possible, 
using the asymptotic expansion of T(A), to determine an infinite set of local, commuting 
Hamiltonians on appropriate Poisson subspaces of £{Q*), or reductions thereof. 

Following the above, let us take r] in (1.7) to be of the form 

T^{X) = 7]o + m\ (1-20) 

and consider the subspace A^o,i C A^o,i consisting of elements of the form 

(;U,e)(A) = (J + CiA, eo + eiA), (1.21) 

where Ci G gin, eo,ei G C are constants and J : S^ ^' gin is arbitrary. Since this is a 
Poisson subspace of i{Q*) with respect to the i?^ Lie-Poisson bracket, the Poisson bracket 
can be restricted to functions that depend only on J, giving 






(1.22) 



The functional differential ^ : 5"^ ^ gin is defined here by the usual formula, 

^MJ + iK)\t=o= [ trK^ \/K:S'^gln. (1.23) 

at J gi oJ 



'SI 
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From now on we choose the values of the Casimirs cq and ei to be cq = 1 and ei = 0. The 
space of Poisson brackets given by (1.22) has a basis given by the choices 77(A) = 1 and by 
rj{X) = A. A common way of describing the above is to say that the space of functions of 
J has two compatible Poisson brackets: 

WM,{J) = - I trCi[g,^] (1.24a) 



and 



(-^«^)-4.M^'^'§'-(^)'^i- <-^^' 



We have denoted the Poisson bracket belonging to 77 = A* as { , }i_|_i in order to be 
consistent with the traditional terminology of KdV systems later. Note that the first 
Poisson bracket is just the Lie derivative of the second Poisson bracket with respect to the 
vector field that generates translations of J in the direction — Ci. 

In general, "interesting" examples result from appropriate further reductions of a space 
M.0^1. A key observation in this respect is that the group consisting of those A-independent 
nonabelian gauge transformations g : S^ —^ Gin for which 

gC^g-' = C^ (1.25) 

is a symmetry group of the AKS system restricted to 7Wo,i- Indeed the transformations 
of 7Vfo,i obtained from (1.17) by using this group preserve both the compatible Poisson 
structures and the monodromy invariants. Thus one can use this group or any of its 
subgroups when searching for "nice symmetry reductions" of the AKS system. 

Remark 1.2. The apparent generalization obtained by considering M._m^n with n + m > 1 
instead of A1o,i does not add any new interesting structure, since the reduction only affects 
Un-i and Un] the other terms remain generic. 

1.2. How to obtain local invariants 

Let us consider the operator C = {d + fx) where fx : S^ ^' £{gln)- There is a fairly 
well-known sufficient condition on the form of the function fu, that one can impose in order 
to guarantee the locality of the monodromy invariants of C This condition involves the 
graded regular elements of the affine algebra i{gln), and we shall explain a variant of it 
below. 

Consider a fixed element A G i{gln) and denote its kernel and image in the adjoint 
representation by 

/C = Ker(adA) and J = Im(adA). (1.26) 
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Of course, /C is a subalgebra of i{gln)- For A a regular element one has 

i{gln) = }C+X, /CnJ={0}, /C : abelian . (1.27) 

We are interested in regular elements that are homogeneous with positive degree with 
respect to some integral grading. An integral grading of £{gln) can be defined by the 
eigenspaces of a linear operator d-N^n '■ ^{gln) -^ ^{gln) of the form 

d 
dN,H = NX—+a.dH, (1.28) 

dX 

where A^ is a non-zero integer and H is a, diagonalizable element of gl^ with integral 
spectrum in the adjoint representation. This formula in fact defines a derivation on i{gln) 
with integral eigenvalues and finite dimensional eigenspaces. 

The following proposition, which generalizes the corresponding well-known result for 
the homogeneous grading case ([1], [19], [22]), states the existence of a solution of the 
linear problem given by a series that can be computed by an algorithm involving only 
linear algebraic operations and integrations. 

Proposition 1.2. Let a grading (1.28) of £{gln) be given. Let A be a regular homogeneous 
element of grade / > 0. Consider a function n : S^ ^' £{gln) of the form 

fi{x) = {q{x) + A) , (1.29a) 

where 

q{x) = J2 q''{x) with dN,H{q'') = kq^ ■ (1.29b) 

k<i 

Then the linear problem 

{d^ + q{x) + A)^x) = (1.30) 

has a unique solution of the form 

^x) = {I + W{x))e^^''\l + W{0))-^^0) , (1.31) 

where 
F{x)elC, W{x)el, and W^(x) = ^ M^^(x) with dN,H{W^) ^ kW^ . (1.32) 

fc<0 

Here the W^ 's are uniquely determined differential polynomials in the components of q and 
F is given by 

F{x) = - [ dy [qK{y) + {q{y)W{y)),c + A] , (1.33) 

Jo 
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where the subscript /C refers to the /C component in the decomposition £{gln) = IC + X. 
Corollary 1.3. The monodromy matrix of C = {d + q + A) is conjugate to 



r.271- 



exp(F(27r)) =exp(- / dx[q!cix) + iq{x)W{x))!c + A]) , (1.34) 

Jo 

and thus its invariants are functions of integrals of local densities in q{x). 

Proof. The procedure is essentially the same as given in a special case in [22]. By substi- 
tuting the ansatz 

^ = {I + W)e^-^, (1.35) 

where / is the unit matrix and \1/ is a constant, into (1.30) we obtain 

W + {I + W)F' + {q + A){I + W) = 0. (1.36) 

If we decompose this equation according to i{gln) = X + /C by using (1.32) together with 
X/C C X and /C/C C /C, then the /C-component gives 

F' + [q{I + W)]K + A = 0, (1.37) 

which (up to a constant) can be integrated to give (1.33), since /C is abelian. By substituting 
(1.37) into (1.36), the X-component gives 

[A,W] + W' - W{q + qW)K: + {q + qW)i = . (1.38) 

One can solve this equation recursively for the VF'^'s by using the grading assumptions of 
(1.29) and (1.32) together with the fact that ad A maps X to X in a one-to-one manner 
since A is regular. This procedure obviously yields the VF'^'s as differential polynomials 
in q. Finally, the integration constant \E' in (1.35) is fixed by the initial condition, giving 

(1.31). D 

Note that the recursive procedure appearing in the proof, combined with the diago- 
nalization of the generators of /C, is also useful for computing the monodromy invariants in 
practice. It should also be noted that in the above we have not considered the convergence 
of the series solution at all. It is well known that such series do not converge in general, 
and are to be considered as asymptotic expansions in A (or alternatively as formal series). 

1.3. The list of graded regular elements 

We have seen that the graded regular elements of i{gln) can be used to impose con- 
straints on the form of C = {d + ^) leading to local monodromy invariants. The sug- 
gestion of De Groot et al [7] was to use the graded regular elements of the inequivalent 
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graded Heisenberg subalgebras of the loop algebras [12] to construct generalizations of the 
Drinfeld-Sokolov hierarchies. The graded Heisenberg subalgebras of the loop algebra £{gln) 
(maximal abelian subalgebras that acquire a central extension in £{gln)^) have been given 
an explicit description recently in [13], where the authors were interested in the related 
vertex operator constructions. By using this description, we shall show that graded regular 
elements exist only in some exceptional Heisenberg subalgebras. The complete list is given 
by Theorem 1.7 at the end of this section. 

The graded Heisenberg subalgebras of £{gln) can be associated to the partitions of n 
in the following way [13]. First, for m any natural number, we define the following m x m 
matrices: 

m — l 

Am = XCm,! + 2_^ ^k,k+l ■> 
fc=l 

TTi — 1 

Hm = diag [i, (j - 1), ... , -{3 - 1), -3] , 3 = -^r- ' (1-39) 



and 



am = exp[2TviHm/m] 



where eij is the standard elementary matrix with entry 1 in the ijth place and elsewhere. 
Let a partition of n be given by 

n — ni + n2 + ■ ■ ■ + rik , where ni > n2 > ■ ■ ■ > n^ > I ■ (1-40) 

We associate to this partition the n x n matrix 

cr = diag [cTn, , (Jna ,•••, cr^ J , (1.41a) 

and denote by N the order of the inner automorphism of gl^ acting through conjugation 
by a. If we let 

TV' = lcm(ni,n2,...,nfc) , (1.41b) 

we have 

f iV', if A^'(^ + ^) is even for all ij: 

1 2N\ if A^'(^ + ^) is odd for some z,j. ^ ' ' 

We then introduce the n x n diagonal matrix H via the equation 

a = exp[27riH/N] , (1.41d) 

and consider the grading of i{gln) given by the eigenvalues of 

d = NX— + 8idH. (1.42) 

aA 
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The Heisenberg subalgebra corresponding to the partition (1.40) is spanned by the n x n 
"block-diagonal" matrices A of the following form: 



A 



y2K 



2 

n2 



VkA^'' 



rik 



(1.43) 



where the li (z = 1, 2, . . . , /c) are arbitrary integers and the i/i are arbitrary numbers. This 
maximal abelian subalgebra of £{gln) is invariant under the grading operator (1.42). An 
element A of the form given by (1.43) is regular if /C = Ker(ad A) is exactly the Heisenberg 
subalgebra (and not a larger space). We next investigate the existence of the graded 
regular elements for some simple partitions, from which we shall then be able to read off 
the answer for the general case. 

The simplest case is that of the homogeneous Heisenberg subalgebra, which belongs 
to the partition n = l + l + - — hi, when d = X-j^ and the graded regular elements are of 
the form 

A = AMiag[yi,y2,...,yn], Vi 7^ Vj . V/c. (1.44) 

The principal Heisenberg subalgebra belongs to the other extreme case when n is "not 
partitioned at all" . In this case d = nX-^ + ad H^ and the graded generators are the powers 
of the "Drinfeld-Sokolov matrix" A„. We have 

^i+mn ^ ^m^i^ ^ rf(A^+""'") = (/ + mn)A^+'"'^ , < / < (n - 1) , Vm . (1.45) 

It is obvious that A^"*""^"^ is regular if and only if A^ is regular. The Drinfeld-Sokolov matrix 
A^ itself is regular since its eigenvalues are the n distinct nth-roots of A, and from this one 
also easily verifies the following by looking at the eigenvalues of A^. 

Lemma 1.4. The element A^ (^1 < / < (n — 1)) is regular if and only if n and I are 
relatively prime. 



Consider now a partition of the type 

n = rii + 712 7 with ni > ?i2 > 1 . 



(1.46) 



Lemma 1.5. In the case (1-46) there is no graded regular element in the Heisenberg 
subalgebra. 

Proof. As candidates for graded regular elements, it is enough to consider the matrices of 
the form 



A 



2/2A: 



n2 



with yiy2 7^ , h ^ mod Ui . 
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(1.47) 



Case (i): Assume that ni and 77-2 are relatively prime. We can check from the definition 
of the grading that there is no graded element of the form (1.47), since 



ni 





'^2 



and 



On, 



A'2 
■n2 



(1.48) 



have different grades. 

Case (ii): If rii and n2 are not relatively prime and m > 1 is their greatest common divisor 

then the graded elements of the form (1.47) are those for which 

fl-i fir) 

h — k — and I2 — k — , k: any integer. (1.49) 

m m 

This implies by Lemma 1.4 that neither of the Aj^. (i = 1, 2) is regular, and therefore there 
is no graded regular element of the type (1.47) either. n 

Let us also note the following rather obvious fact. 
Lemma 1.6. There is no graded regular element in the Heisenberg subalgebra if the 
partition is of the type 

n = m + l-i hi, l<m<(?i-l), (1.50) 

i.e. if it consists of a "non-singlet", m, and more than one "singlets" (I's). 

It follows from the "block structure" of the Heisenberg subalgebras given by (1.39)- 
(1.43) above that graded regular elements do not exist for any of those partitions which 
contain a subset of the type appearing in Lemmas 1.5 and 1.6. Hence the only cases not 
excluded are the "partitions into equal blocks", 

p times 

77, = pr = r + --- + r, (1-51) 

and the cases "equal blocks plus a singlet" , 

p times 

n = pr + l = r + --- + r+l. (1-52) 

On the other hand, by inspecting the eigenvalues of the generators of the Heisenberg 
subalgebras, we can establish that graded regular elements do indeed exist in the cases 
(1.51), (1.52), and the following theorem gives the complete list. 
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Theorem 1.7. Graded regular elements exist only in those Heisenherg subalgebras of£{gln) 
which belong to the special partitions (1.51) or (1.52). In the equal block case (1.51) with 
r > 1 the graded regular elements are of the form 

2/2 Aj^ 



A = A^ 



ypK[. 



(1.53) 



where 

l</<(r-l), ?/i^0, yl^y] i,j = l,...,p, i^j, 

with I relatively prime to r and m any integer. The element A is of grade (/ + mr) with 
respect to the grading operator given by (1.4^) where N = r and 



p times 



H = diag[Hr,Hr,...,Hr]. (1.54) 

In the equal-blocks-plus-singlet case (1.52), the graded regular elements are those n x n 
matrices which contain an (n — 1) x (n — 1) block of the form given by (1.53) in the 
"top-left corner" and an arbitrary entry in the "lower-right corner". The relevant grading 
operator is given by (1.42) with N = r, 

p times 

H = diaig[IIr, Hr, . . . , Hr^O] if r is odd] (1.55a) 

and with N = 2r, 

p times 
. ^ . 



H = diag[2Hr,2Hr,...,2Hr,0] ifriseven. (1.55b) 

Remark 1.3. Let us designate the ordered eigenvector basis of the matrix H in (1.54) as 

"^i,l'"^J-l,l' • • • 7-^-j,i;-^j,27-^j-l,2, • • • j^-j,2] ■ ■ ■ '■>^j,p^^j-l,p^ ■ ■ ■ ^^-j,p- (1.56) 

Here j = (r — l)/2, the first index is the eigenvalue and the second one orders the r x r 
blocks. It is often convenient to use the re-ordered basis 

"^i,i'"^j,25 • • • 7-^j,p;-^i-i,i'^j-i,2, • • • ,Xj-i,p', ■ ■ ■ ',^-j,i,^-j,2, ■ ■ ■ ^^-j,p- (1-57) 

When expressed in the new basis the matrix H of (1.54) may be written as Hj, ^ Ip and 
the graded regular element A given by (1.53) takes the form 

X'^Al(»D, where D ^ diag{yi,y2, . . . ,yp) . (1.58) 

In the following section we shall consider reductions based on the grade 1 regular 
elements belonging to the equal block case. 
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2. Equal block reduction of the AKS system 
to the matrix Gelfand-Dickey hierarchy 

We have seen that graded regular elements exist only in those Heisenberg subalgebras 
of£{gln) which correspond to partitions into equal blocks or equal blocks plus a singlet. The 
purpose of this section is to study in some detail symmetry reductions of the general AKS 
system that are based upon the grade 1 regular elements of the Heisenberg subalgebras 
belonging to the partitions into equal blocks, n = pr, generalizing the p = 1 case considered 
in [1]. The final result of the analysis below is that the reduction of the bihamiltonian 
manifold Ai := A^o,i carrying the commuting family of ad*-invariant Hamiltonians yields 
the pxp matrix version of the well-known Gelfand-Dickey r-KdV hierarchy. More exactly, 
we shall establish the following: 

1. The reduced phase space is the space of rth order, pxp matrix differential operators 
carrying the first and second Gelfand-Dickey Poisson brackets. The second Poisson 
bracket algebra is an example of a classical W-algebra. 

2. The commuting hierarchy of Hamiltonians resulting from the monodromy invariants 
is given by the componentwise residues of the fractional (including integral) powers 
of the pseudo- differential operators obtained by diagonalizing the matrix differential 
operators. 

The exact statements are given by Theorem 2.4 and Corollary 2.11 in Sections 2.2 and 
2.3. These results generalize the analogous results proven by Drinfeld and Sokolov for the 
scalar case p = 1. We shall in fact use many of their methods, but at the same time 
introduce some simplifications (at least to our taste) in the proofs. 

Remark 2.1. For p > 1 the subhierarchy provided by the trace- residues of the fractional 
powers of the matrix differential operators is not exhaustive, since it does not include 
the Hamiltonians obtained from the integral powers of the corresponding diagonal pseudo- 
differential operators, which also appear in Corollary 2.11. 

2.1. A local symmetry reduction of the AKS system 

After reordering the basis as explained previously, the generators of our Heisenberg 
subalgebra are the n x n matrices of the form A^ (E) D, where Aj; is the k^^ power of the 
r X r DS matrix, k is an arbitrary integer, and D is an arbitrary pxp diagonal matrix. 
The generator A^ eg) D is of grade k under the grading defined hy d = rX-^ + ad H with 

H = Hr^lp = diag [jlp, (j - l)lp, . .. , -(j - l)lp, -jlp] , j = ^^^. (2.1) 
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Choose a grade 1 regular element oi the Heisenberg algebra, (cf. (1.53), (1.58)); i.e., an 
element 

K:=Kr®D (2.2) 

such that D^ has distinct, non-zero eigenvalues. Define the A-independent constant ma- 
trices Co and Ci by the equality 



K = Kr®D •=Cq + ACi 



(2.3a) 



More explicitly, these matrices are given in block form by 



A 






D 


■ 


■ 







D 


. 









. D 


kD 







■ 



Co + ACi 



(2.3b) 



We start reducing the general AKS system carried by l{gln )* by confining ourselves to 
the space M. consisting of elements of the form 



£ = a + J + ACi 



J : S^ ^ gin . 



(2.4) 



This is an example of a Poisson subspace of the type A^o,i considered in Sec. 1.1, and 
therefore it carries the two compatible Poisson brackets given by (1.24a,b). 
Consider the decomposition 



g^n = gin + gin + gin 



+ 



(2.5) 



induced by the eigenvalues of ad if, eq. (2.1), where the summands are the subalgebras 
of block lower triangular, block diagonal, and block upper triangular matrices (with p x p 
blocks), respectively. On account of the relation 



[gl- , Ci] = 



the group A/" of transformations 



^J -.C^ e^Ce-f , with / : 5^ -> gl; 



(2.6) 



(2.7) 



is a symmetry group of the AKS system carried by M.; i.e., these transformations preserve 
the two Poisson structures and the monodromy invariants. Next we define a symmetry 
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reduction of the AKS system using Af in such a way as to ensure the locality of the reduced 
system. That is, consider the following two step reduction process, which is an obvious 
generalization of the one used in [1]. First, we restrict our system to the "constrained 
manifold" Aic C Ai, defined as the set of £'s of the following special form: 

C = d+{q + Co)+XCi = d + q + A, q:S^^{gl-+gll) . (2.8) 

Here Co is the constant matrix given by (2.3), and q is required to be block lower triangular. 
(Note that with respect to the second Poisson structure (1.24b), this is just fixing a level 
set of the moment map generating the Hamiltonian group action (2.7), and the reduction 
procedure is essentially that of Marsden-Weinstein at the level of a Poisson manifold rather 
than a symplectic one. With respect to the first Poisson structure (1.24a), the constrained 
quantities determining the form of £ in (2.8) are all Casimirs and hence this just determines 
a Poisson submanifold.) Second, we factorize this constrained manifold by the symmetry 
group A/", defining the reduced phase space 

MreA = Mc/M . (2.9) 

To put it another way, we factorize out the "gauge transformations" generated by N by 
declaring that only the A/'-invariant functions of C are physical. The nice features of this 
reduction are that 

i) The monodromy invariants oi C E M.c can be computed algebraically as asymptotic 
series in A which depend on q through integrals of local densities formed from its 
components and their derivatives. 

ii) The compatible Poisson structures on M. induce compatible Poisson structures on 

MreA- 

iii) The gauge orbits in M.^ allow for global, differential polynomial gauge sections, which 
give rise to complete sets of gauge invariant differential polynomials. 

Statement i) follows immediately since we have chosen A^c so that the conditions of 
Proposition 1.2 are satisfied. Statement ii) means that the compatible Poisson brackets 
carried by M. can be consistently restricted to the gauge invariant functions in C°°(A4c), 
whose space can be naturally identified with C°°(A/1 red)- This can be seen as a consequence 
of the Dirac theory of reduction by constraints as follows. We first note that, by choosing 
some basis {7^} of gl~^ the constraints defining M.^ (Z AA can be written as 

Xi{x) = where Xi{x) = tr7i(J(x) - Co) . (2.10a) 
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It is easy to verify that these constraints are first class; i.e., 

{x^{x) . Xk{y)}\M. = , (2.10b) 

for any of the compatible Poisson brackets on Ai. We next notice that the functions Xi{^) 
are the generating densities (i.e., components of the moment map) of the A/" symmetry 
transformations with respect to the second Poisson bracket (1.24b). Thus the theory of re- 
duction by constraints (which in this case is just the Poisson version of Marsden-Weinstein 
reduction) teUs us to factorize the constrained manifold by these transformations; the sec- 
ond Poisson bracket algebra closes on the gauge invariant functions on A^c, inducing a 
Poisson structure on the factor space. On the other hand, the Xi do not generate any 
transformations on A1 under the first Poisson bracket (1.24a); i.e., they are Casimir func- 
tions. Therefore the first Poisson bracket can in principle already be restricted to C°°{Aic) 
without any factorization by A/". Then A/" becomes a group of Poisson maps with respect to 
the restricted bracket, which can further be reduced to a Poisson bracket on the invariant 
functions. In this way, we naturally obtain two induced Poisson brackets on A/1 red from 
those on Ai , and the induced Poisson brackets are compatible because the original brackets 
(1.24a,b) were compatible. 

One is always interested in gauge invariant objects and convenient gauge fixings when 
describing systems with gauge symmetries. In the present example, as in the p — 1 case 
of [1], a convenient gauge section is defined by the subspace 

r 
V -.^{d + Y^ er,^ ® Vr-i+l + A \ Vk : S^ ^ glp} , (2.11) 

i=l 

i.e., the space of £'s in which the matrix function q is allowed to have non-zero {p x p 
block) entries only in the last row. This global section of the gauge orbits can be reached 
from an arbitrary point C = (d + q + A) gA/1c by a unique gauge transformation that 
depends on g in a differential polynomial way. It follows that the components of the Vi 
provide a basis (generating set) for the gauge invariant differential polynomials which can 
be formed from the components of q. 

Remark 2.2. A distinguished gauge invariant differential polynomial is obtained by re- 
stricting 

VH = \tr{J^)+tr{HJ') (2.12) 

from ]v[ to A^c- The density Vh satisfies the Virasoro algebra under the second Poisson 
bracket, and therefore it generates an action of Diff(5'^) on M. that survives the present 
reduction. Since it contains this Virasoro density, the second Poisson bracket algebra of the 
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gauge invariant differential polynomials can be regarded as an extended conformal algebra; 
that is, a classical W-algebra. This W-algebra is a member of a natural family of extended 
conformal algebras that recently has been studied in [4], [5]. The reader is referred to [5] 
for a detailed description of differential polynomial gauge fixings, like the gauge (2.11), 
and for the related construction of a generating set for the W-algebra consisting of Vh and 
conformal tensors. 

In the next section we shall also need the "block diagonal gauge" given by 

r 

e:={d + Y^e^,^^9, + A\9,:S^^glp}. (2.13) 

It should be noted that G defines only a partial gauge fixing, which has a finite dimensional 
residual gauge freedom. Of course it is nevertheless true that any gauge invariant element 
of C'^{Aic) can be recovered from its restriction to G. In particular, the Poisson bracket 
of any two gauge invariant functions in C°° (Aic) can be recovered from its restriction to G, 
and the nice feature is that for the second Poisson bracket this restriction can be computed 
in terms of the "free current algebra" of the ^i's. For later reference we summarize this 
fact as a lemma. The proof is similar to that for the p = 1 case [1] . 

Lemma 2.1. Let ip^ifj E C°°{M.c) be gauge invariant functions and consider ^ = {</?, 'i/j}2 G 
C°^{Aic), which is well defined and is also gauge invariant. Let us denote the restriction 
of these functions to G by if>,^ and ^, respectively. Then we have 



We consider the decomposition J = J_ + Jq + J-|- defined by means of (2.5) and write 
out formula (1.24b) in terms of the partial functional derivatives corresponding to these 
variables. We then compute the value of {</?, 'i/'}2 at an arbitrary point on the "gauge 
slice" G of the block diagonal gauge, where J_ = 0, Jq = diag [6*1, . . . , 6*^,], J_|_ = Cq. 
(For this computation we can extend </? and t(j from Aic io Ai in an arbitrary way since 
they are invariant under the gauge transformations generated by the first class constraints 
specifying this constrained manifold.) Using the grading structure, we find that 

\ (2.15) 



i=l -^ ^ 
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where ^, -^ are understood as the p x p blocks constituting the block diagonal jj-, jj-. 
This immediately gives (2.14) since at the point J = Jq + Cq we have ^ = ^, and 
similarly for t/j. n 

Now let M be the space oi p x p "matrix Lax operators" of the following form: 

r 

L = {-D-^Yd'' + ^u,d''-\ ur.S^^glp. (2.16) 

As in the p = 1 case [1], M can serve as a natural model of the reduced space M.c/N'. To 
see this let us consider the linear problem 



£* = £ 



'i/'2 



0, (2.17) 



for C given by (2.8), where the entries il^i of \1/ are p-component column vectors. This 
system of equations is covariant under (2.7) if we let \E' transform as 

^ -^e^* . (2.18) 

Denoting the zjth p x p block of q by qij, the system (2.17) can be recast in the 
following form: 

V'2=(-I^)-'(a + Qll)V'l 



where L is an operator of the form (2.16) whose coefficients ui are uniquely determined 
differential polynomials in q. Notice now that the component V'l is invariant under (2.18), 
because / is a strictly block lower triangular matrix. This implies that the potentials 
Ui = Ui[q] entering in Li/ji = Xi/ji must also be invariant functions of q under (2.18). In 
other words, the operator L attached to C by the equivalence of (2.17) and (2.19) actually 
depends only on the AT-orbit of C in A4c- Thus it gives rise to a map m : A4red ^ M . It 
is easy to see that this is a one-to-one map. The inverse m~^ can be given by attaching 
to an arbitrary operator L in (2.16) the unique orbit in Aic which intersects the gauge 
section V (2.11) in the point 

r 

C = d + ^er,t^Vr-r+i+A where Vi = {-ly+^-'Du^D"^-' . (2.20) 
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This equation provides an identification of tfie space M witli tfie gauge section V and 
fience witfi tfie reduced pliase space A^red, 

M ~ y ~ Mred ■ (2.21) 

It follows from the above that the natural mapping tt : G ^ M is given by the 
factorization formula: 

L = {-D-^fd' + uid'-^ ^■■■ + Ur = {-D-\d + 9r)) ■ ■ ■ {-D-\d + ^i)) , (2.22) 

generalizing the scalar case ([23], [1]). Lemma 2.1 can be reformulated as saying that 
TT, which is the usual Miura map when p = 1, is a Poisson mapping if the space G 
carries the Poisson bracket appearing on the right hand side of (2.14) and M carries the 
Poisson bracket induced from the second Poisson bracket on Ai via the reduction and the 
identification M ~ TWred. 

2.2. The Gelfand-Dickey form of the reduced Poisson structures 

In this section we prove a theorem that establishes the equivalence of the Poisson 
structures naturally carried by the spaces TWred and M. This theorem will rely on the 
preliminary Lemma 2.2 which follows*. This concerns the Poisson Lie group property of 
the so-called Sklyanin bracket as pointed out by Semenov-Tian-Shansky [20]. 

Let A be an associative algebra. Suppose that Tr : ^ ^ C is a non-degenerate 
trace-form on A; i.e., a linear mapping with the property that the formula (a, b) = Ttab 
defines a non-degenerate, symmetric bilinear form on A. As usual, identify the space A 
with its dual A* (or a subspace thereof) by means of the pairing defined by Tr . Suppose 
that A and B are disjoint subalgebras of A such that A = A + B and that with respect to 
the pairing ( , ) we have A-^ = A and B-^ = B. Let Pa and Pb be the projection maps 
on A defined by the splitting A = A + B. Since A and B turn into Lie subalgebras of 
A with respect to the natural Lie algebra structure on A given by [a, b] = ab — ba, then 
R = Pa — Pb is a skew-symmetric r-matrix on A, and the AKS construction applies here 
too. 

For any function </? G C°°{A), we define Va</7, its gradient at the point a G ^, by 

^V'(a + t6)|t=o = (&,V„v?) = Tr6VaV?, \/beA. (2.23) 

In addition to the R Lie-Poisson bracket given by 

W,ilj}^^Ha) := Tr {a[X,Y]R) , where X = V^v?, ^ = V.V' , (2.24) 



This lemma was explained to us by A. G. Reyman. 
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we can also define a second Poisson bracket on A: 

{(/?, V}^^^(a) := ^Tr {YaR{Xa) - aYR{aX)) . (2.25) 

This Poisson bracket is known as the quadratic r-bracket, the Sklyanin bracket, or the 
second Gelfand-Dickey bracket. The proof of the Jacobi identity for the Poisson bracket 
(2.25) is easy if we remember that because of the derivation property it is sufficient to 
check it for Hnear functions; i.e., functions of the form /^(a) = Tr (aX), X ^ A. In the 
proof we also use the fact that R satisfies the modified Yang-Baxter equation, 

[RX,RY]=R{[RX,Y] + [X,RY])-[X,Y] yX,YeA. (2.26) 

Note that the R Lie-Poisson bracket and the quadratic r-bracket are compatible. 

Lemma 2.2. With respect to the quadratic r-bracket on A, the mapping m : A x A ^' A 
given by m{a, b) = ab is a Poisson mapping when Ax A has the product Poisson structure 
corresponding to the quadratic r-bracket on each of its components. 

Proof. For a pair of elements X,Y E A, consider the linear functions on A given by 

/,(/)= Tr/X, /^(/)=Tr/y. 

Then 

m* f^{a,b) =TrabX and m* fy. {a, b) = Tr abY . 

Variation of m* f^ with respect to a for fixed b gives 

VaTO*/xl(a,6) = bX , 

similarly 

V6m*/xl(a,&) = Xa . 

Thus 



{m*/x,^*/y}!4x^(a,^) =oTr {bYaR{bXa) - abYR{abX)) 



+ ^Tr {YabR{Xab) - bYaR{bXa)) 



:-Tr {YabR{Xab)-abYR{abX)) = (m*{/^, /^jjl^) (a, 6) 
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We shall apply this lemma to the following example: 

N 

A={X= J2 ^sd'l^Xs-. S^ ^pxp matrices, \fN>0}:=GDM, (2.27) 
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the space of pseudo-differential operators with px p matrix coefficients. We caU this space 
GDm for "matrix Gelfand-Dickey" . Multiphcation of matrix pseudo-differential operators 
is defined in the usual way ([14], [16]). The trace form, Tr : GDm -^ C, is given by 



TrX:=/ trres(X)=/ tr X_i , (2.28) 

where tr is the ordinary matrix trace. The two subalgebras A and B in this case are given 
by 

N -1 

A:^GDM^={X = Y,Xsd'}, B:^GDm_={X= J^ X,c>^ } (2.29) 

s=0 s= — oo 

It is usual to write Pa{X) = X+ and Pb{X) = X-. One sees by inspection that for any 
/c > 0, the space 

Mo,fc :={X = XqO^ + Xid^-^ + ■ ■ ■ + Xfc I Xo fixed } C GDm (2.30) 

is a Poisson subspace with respect to both Poisson brackets (2.24) and (2.25). 

In particular, the space M of operators L defined by (2.16) is such a subspace, on 
which the two Poisson brackets take the form 

{V?,^}(1)(L) = / trres (L[y_,X_]) , (2.31) 

{^,i)}^'^\L) = I trres {YL{XL)+ - LY{LX)+) , (2.32) 

where X := Vl(/7, Y := V^^ for Lp, i/j E C°^{M). It is now easy to prove the following 
corollary of Lemma 2.2. 

Corollary 2.3. The mapping tt : ^ M given by eq. (2.22) is a Poisson mapping where 
the Poisson bracket of f^h E C°°(0) is given by 

,/..,(.......,^|:/^,.(«4,|,.(|)'|). (.33, 

and the Poisson bracket of (pji/^ G C°° (M) is given by the second Gelfand-Dickey bracket 
(2.32) on the manifold M. 

Proof. It is sufficient to check that on the Poisson subspace {—D~^{d + 6)} of GDm (a 
subspace of type Mo,i) the Poisson bracket (2.25) becomes just the r = 1 case of (2.33), 
that is, if v?, ^ G C^{{e}), then 
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We may now state the main result of this section. 

Theorem 2.4. Under the identification A^red — M given in (2.21), the Poisson brackets 
on the space TWred = M-c/M obtained by reduction from the Poisson brackets (1.24a) and 
(1.24b) on M. are equal respectively to the first and second G elf and- Dickey Poisson brackets 
given by (2.31) and (2.32). 

Proof. Consider the one parameter group of transformations on M. defined by 

g^:C^{C- tCi) , r e R (2.35) 

These transformations preserve Aic '^ M. and commute with the action of A/". Thus 
we have a corresponding one parameter group of transformations {^r} on M ~ Aic/-^^ 
operating as 

Qr : {L \-^ L + rip) (Ip : the p x p unit matrix) . (2.36) 

We now recaU that { , }i given by (1.24a) is the Lie derivative of { , }2 given by (1.24b) 
with respect to the vector field generating the fiow (2.35). It follows that a similar relation 
holds for the corresponding induced brackets on M ~ Aired with respect to the generator 
of the projected fiow (2.36). On the other hand, one also sees by inspection that the 
first Gelfand-Dickey bracket on M; i.e., the restriction of (2.24) given by (2.31), is the 
Lie derivative of the second Gelfand-Dickey bracket (2.32) with respect to the generator 
of (2.36). Therefore the equality of the "first brackets" claimed in the theorem follows if 
we prove the equality of the "second brackets" . This latter is obtained as an immediate 
consequence of Lemma 2.1 of and the above corollary to Lemma 2.2. n 

Remark 2.3. It is clear from the proof that Theorem 2.4 remains valid if we let D be any 
invertible p x p matrix in the construction. The choice of a graded regular element for 
A in (2.2) is made in order to guarantee the existence of local monodromy invariants by 
Proposition 1.2. In the next section this assumption will be fully used in determining the 
Hamiltonians of the reduced AKS system. 

Remark 2.4. Let Re : GDm — * GDm be the mapping defined by right-multiplication by 
the invertible (constant) matrix c; i.e., Rc{ci) = ac for any a G GDm- This is a one- 
to-one mapping that preserves the second Gelfand-Dickey bracket and carries the first 
Gelfand-Dickey bracket into the Poisson bracket { , }c on GDm given by 

{ip,ij}'^^\a) = TT{ac-^[cX,cY]R), where X = Va^,Y = Va^, (2.37) 

which is compatible with { , }^'^> for any c. By means of this mapping, with c = {—DY, 
we can map the space M of L's in eq. (2.16) onto the space M consisting of the rth order 

26 



differential operators of the form 



1 = 8"^ + ^Uid' 



Ui-. S^ ^ gl 



p ■ 



(2.38) 



This is the reduced space we would have obtained by replacing the matrix D by the 
unit matrix in the reduction procedure described in the previous section. (The results 
of Theorem 2.10 and Corollary 2.11 below, which characterize the reduced monodromy 
invariants, would not then be applicable.) The one-to-one mapping between M and M 
provided by Re can be made into an isomorphism of bihamiltonian manifolds if one lets 
M C GDm carry the restrictions of the two Gelfand-Dickey brackets, while M C GDm 



carries the restriction of the "modified first bracket" { 
restriction of the second Gelfand-Dickey bracket. 



}^P with c = {-Df, plus the 



2.3. Computation of the Hamiltonians 

In the previous section, we identified the Poisson brackets carried by the reduced space 
A^i-ed — M as the first and second Gelfand-Dickey Poisson brackets. In this section we 
establish a description of the Hamiltonians of the reduced AKS system in terms of the 
reduced-space variables -ui, . . . , Ur- This will result from a computation of the eigenvalues 
of the monodromy matrix of the linear problem £\1/ = for an arbitrary C &V; i.e., for 

/O ■■■ 0\ 

£ = a + Co + ACi + 







\Vr 







Define the matrix A by 



A := {-D) 



(2.39) 



(2.40a) 



Note that A^ is a nondegenerate and invertible diagonal matrix, by eq. (2.2). Let 



L = A^a'^ + wi^'^-^ + • ■ ■ + w 



(2.40b) 



be the element of M corresponding to C E V . If {^a}a=i,...,rp is a complete set of 
independent (p-component column vector) solutions to 



L(j) = X(f) . 



(2.41) 



define the n x n matrix $ by 



$ 



91 



VA'^-Vl'"'^ 



''rp 



M'r 



rp 



-W(r-l) 
<^rp 



(2.42) 
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Then the columns of $ are a complete set of solutions to £\1/ = and T = $(27r)$(0)~^ 
is the monodromy matrix, whose invariants will be computed. 

Since A*^ is nondegenerate and diagonal, we can, through the usual, recursive ap- 
proach, find an operator g of the form 



CXD 



g = I + J29^^-\ (2.43) 

i=l 

with gi{x + 2%) = gi{x) for all i, such that 

L = gtg-^ (2.44) 

with L a diagonal matrix pseudo-differential operator; i.e., with L of the form 

CXD 

L = A'^d^ + 22 <^i9^~^ 7 <^i '■ ^11 diagonal matrices. (2.45) 

For example, if we require the (/j's to be off-diagonal then we can recursively uniquely 
determine both the (/j's and the a^'s as differential polynomials in the w^'s, by comparing 
the two sides of Lg = gL term-by-term according to powers of d. Later we shall need that 

ai = [ui]diag- (2.46) 

Fix —( to be any rth root of A, 

(-C)'^ = A. (2.47) 

Consider the p x p diagonal matrix asymptotic series 

Vi(x, C)=e'^^"(^(x) + C"'<^i(^) + ---) for C~oo (2.48) 

satisfying the equation 

{U)ixX) = i-CYH^X) (2.49) 

in the following sense (see e.g. [24]). Using the definition 

(a-^e*^^") := (CI^)-"e'^^" (2.50) 

for s any integer, and extending this in the obvious way to pseudo-differential operators, 
we can write 

V'(x, = {Ve'^^'n , (2.51a) 
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where 

V = do{x) + di{x)d-^ + ■ ■ ■ (2.51b) 

is a uniquely determined diagonal pseudo-differential operator. Then the action of any 
pseudo-differential operator P G GDm on ^ip is defined by multiplication of pseudo- 
differential operators, 

(PV')(^, C):=((i'^)e^^"). (2.52) 

The left hand side of equation (2.49) is understood in this sense. If we assume that 
det(/)o(x) 7^ in (2.48), then ip is uniquely determined up to multiplication by an x- 
independent diagonal matrix of the form 

c(() = co+C"^ciH , with detco^O. (2.53) 

At this point we make use of a result given in [1] (Theorem 2.9, pages 1986-7). The 
proposition below states the result in terms of the diagonal pseudo-differential operator L, 
and therefore we are dealing with p different copies of the scalar case. 

Proposition 2.5. If L is of the form (2.45), d may be expressed as: 

oo 

d = -DL\'' + Y, F,{-DL\''')-' , {-D = A-i) , (2.54) 

with 

r.27r 



1 I / \ 

Fo = — {-Dfai and j (kFk + {-D)''res{L'l^'^)j=0, for k > 0, 

'^ 1 /r 

where L^' is the unique diagonal pseudo- differential operator such that 



(2.55) 



L;/'^ = Aa + ^6,a-^ and (L^'^)^ = L . (2.56) 



A further argument from [1] can be applied directly to deduce from (2.49) 

Lemma 2.6. 

(L;/'^V')(^,C) = -C^(^,C)- (2.57) 

Applying d io ip and making use of Proposition 2.5 and Lemma 2.6, we get 

^'= Ud + Y.F,{CD)A^. (2.58) 
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This equation can be solved, giving 

(Dx + J^iCD)-" J^ FkU{0,O- (2.59) 

It follows that 

iP{x + 2nX) = ^ixXh{C), (2.60) 

where 

2nCD + J2iCDr' J^ FkjmO- (2.61) 

Define the p x p matrix asymptotic series t/j^x^Q by 

iPixX) := {9i^){xX) , (2.62) 

where g is the pseudo-differential operator appearing in (2.43). We then have the self- 
evident 

Lemma 2.7. The columns of ifj defined by (2.62) are solutions of (2.41)- 

If we expand g in descending powers of —L^ instead of descending powers of d we 
have 

^ = I + mi(-L;/'^)-i + m2(-L;/0"' + ■ • • , (2.63) 

for some rMi,rM2, . . ., with mi{x + 27r) = mi{x) for all i. Then (2.57) gives 

'0(a^, C) = (fi'V') {x, C) = m{x, O^ix, C) , (2.64a) 

where 

m{x, C) := (/ + mi{x)C~^ + m2{x)C~^ -\ ). (2.64b) 

From this we obtain 

Lemma 2.8. 

V'(x + 27r,C) = ^(a;,C)7(C). (2.65) 

Now consider the set of r independent roots of A, given by 

{-a = -e'"^'/'^C}fc=o,...,r-i, (2.66) 

and define ifji by 

ipi{x,C) ■=i^{xXi) ='m{xXi)^{X:Ci) z = 0, l,...,r- 1. (2.67) 
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The point of this construction is: 

Lemma 2.9. 

The columns of the n x n matrix 



$ 



/ V'o 



vA-vr^^ 



V'r-l \ 

A<_i 



A'^-VtV^ / 



(2.68) 



are a complete set of solutions to £\1/ = 0. 



Proof. We only have to prove that det $ 7^ 0. This follows by checking the leading term 
for C ~ 00. Using (2.48), (2.64a,b) and (2.67) we have 



$(x,C) 



/ 1 

-Co 



V(-Co) 



r-l 



1 



\ 



,r-l 



(-C.-i)^-V 



fe'^'^^'^Mx) 



V 



^Cr-lDx 



\ 



(poix)/ 



(2.69) 



for C ~ 00. The determinant of the matrix on the right hand side of equation (2.69) is 
non-zero. n 



From Lemma 2.8 we obtain 



$(a: + 27r,C) = $(a:,C)r(C), 



(2.70) 



where 



r(C)=diag(7(Co),...,7(Cr-i)). (2.71) 

It follows that T = $(27r)$(0)~^ is conjugate to F. By combining this with (2.61), (2.55) 
and (2.46), we arrive at the main result of this section. 

Theorem 2.10. The monodromy matrix of C is conjugate to the diagonal matrix T given 
by 



/27iCoD-ln{Co) 



T = exp 



where 



2nC^D - mC^) 



\ 



2T:Cr-lD-ln{Cr-l)j 



(2.72) 



/2-K 00 „27r 

{-DY[uiW, + Y.^-C)-\ res(Lt/^ 

fc=i •^'^ 



(2.73) 
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and L is the diagonalized form of the operator L E M corresponding to C E M.c- 

Corollary 2.11. All Hamiltonians of the reduced AKS system carried by M are generated 
by the ones in the following list: 

(2.74) 
The corresponding flows are subject to the relations 

{f,nkA^^^ = {f,nk+rA^^^ ^fec^{M), yi,k. (2.75) 

The first Hamiltonians in each of the bihamiltonian ladders belonging to fixed k mod r and 
fixed i are Casimirs of the first G elf and- Dickey bracket, 

{/,^m}^'^=0, V/gC°°(M) for fc = 0,l,...,(r-l). (2.76) 

The number of independent bihamiltonian ladders in (2.75) is pr — 1 = n — 1 since 

p 

y^ T-Cmr,! = , for any m = 1, 2, . . . , (2.77) 



i=l 



and X]r=i '^0,1 is a Casimir with respect to both G elf and- Dickey Poisson brackets. 

Proof of Corollary 2.11. It is clear from the form of the matrix T in (2.72) that the 
Hamiltonians of the reduced AKS system; i.e., those obtained by reduction from the ad*- 
invariant Hamiltonians on the dual of the loop algebra, are generated by those in the list 
(2.74). The relations given by (2.75) and (2.76) can be traced back to the "general recursion 
relation" given by (1.10), but can also be verified directly by using the expressions of the 
Gelfand-Dickey brackets, (2.31), (2.32), and the Hamiltonians (2.74). That X]i'=i ^o,i is 
a Casimir of both Poisson brackets also follows by direct verification. As for (2.77), note 
that 

P p2-K j-2-K 

y"Hmr,*=/ trresL"^=/ trresL"^ = 0, (2.78) 

i^i ' Jo Jo 

where the second equality follows from (2.44) and the third equality holds because L is a 
differential operator. n 

Remark 2.5. Let p he a, p x p diagonal matrix whose entries are all rth roots of 1, 

p = diag(pi,...,pp) with {p^y = l. (2.79) 
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There is a unique rth root of L whose leading term is given by pAd. We denote this 
pseudo-differential operator by Lp . By using (2.43), (2.44) and (2.56) we can write 

Lj/'^ = gplY'-g-' , (2.80) 

which implies that 



iXr {L^n = ^Tr {p'^L',/^) = J^ p'^k, • (2.81) 

i=l 

For the Hamiltonian 7i^ := ^TrLp (k ^ mod r) we have 

VlK = ^^'"~^ ■ (2-82) 

From this we obtain, as in the scalar case [1], 

^ = {^ , K}^'^ = {^ , K+J^'^ = [(^p/'^)+ , ^] • (2-83) 

There are r^ possible choices of the matrix p in (2.79) and it is possible to single out p of 
them in such a way that the corresponding Ti^'s form a basis for the linear space spanned 
by the ?ifc,i's, for any fixed /c 7^ mod r. Thus for /c 7^ mod r the relations (2.75) and 
(2.76) follow from (2.83). 

Remark 2.6. Note that the diagonal terms of the potential ui of L in (2.40b) are constant 
along the flows of the hierarchy. This is obvious for flows of the type (2.83), and follows for 
all the flows generated by the ad*-invariants by combining the following two facts, which 
are easy to verify. First, the function 

HKiL) = i-iy / tr [KDuiD^-^) , (2.84) 

where K is any glp valued function on 5'^, generates via the second Gelfand-Dickey bracket 
the following one parameter group of transformations: 

ft-.L^ e^^Le-^^ . (2.85) 

Second, these transformations leave the Hamiltonians (2.74) invariant for arbitrary diago- 
nal K. For example, ii K — diag (i^i, . . . , Kp) is a constant diagonal matrix then 

p 
nK = Y.K,no,^ (2.86) 

1=1 
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is one of the Hamiltonians of the hierarchy. It generates the evolution equation 

L = {L,nKV^^ = [K,L], (2.87) 

whose flow (2.85) indeed leaves [wijdiag invariant (but not the full matrix ui). Finally, we 
could also consider symmetry reductions of the hierarchy under the abelian group action 
given by (2.85) with diagonal K^s. The simplest reduction would be defined by setting 
tr (D^ui) = 0. In this case, the reduced system is the same as the one obtained by using the 
Lie algebra sin instead of gin throughout the construction. Another symmetry reduction, 
which in a sense is maximal, is obtained by setting [•uijaiag = 0. This leads to the system 
studied previously by Gelfand-Dickey [15], Manin [16] and Wilson [17]. 

2.4. Some explicit formulae in the case r = 2 

Here we consider the simplest example r = 2 in order to make our general results 
more concrete. We shall display the explicit form of the Poisson brackets and the first two 
Hamiltonians for each bihamiltonian ladder in our reduced AKS system. 

Denote the general element of the phase space M as 

L = A'^d'^ + ud + w. (2.88) 

We can easily write out the evolution equations generated by an arbitrary 7i G C'^{M) 
via either of the two Gelfand-Dickey Poisson brackets by using 

Vz.H = a-2^ + a-i^. (2.89) 

ou ow 

The Hamiltonian equation generated by means of the first Poisson bracket (2.31) is 
given by 

ti = {t,,7^}(i) = [A2, ^] (2.90a) 



dw 



and 



w 



(-«''"-^^(^)'-(S)'^^-[^^'fi-i-^i' (-«''' 



Observe that [w]diag is in the centre of the first Poisson bracket since it does not change 
under any Hamiltonian flow. In the scalar case p = 1 all the commutator terms drop out 
of (2.90), and for ti = we would recover the first Poisson structure of the standard KdV 
hierarchy. 
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The Hamiltonian equation generated by Ti through the second Poisson bracket (2.32) 
is given by 

\ ou ou I \ ow ow I 

/ / // (2.91a) 

OW J \ou J \ow 



and 



w 






+ I A -^w — tw-^A + -u- — w — w- — u ] — u { -^u 
ou ou J \ ow ow J \ ow 



+H^J ^-^W) -"U'^" 



(2.91b) 



By straightforward computation from (2.43-45) we obtain 

/'27r 1*211 

U2,i{u,w)= / res(L)ii= / /i2,j , (2.92a) 

Jo Jo 



with 



_ >r^ Af + A| , >r^ UikUkiUii 

"-2,^ - 2^ /AS _ /\2\2^ikUki - 2_^ 



^^.(AJ-Ai)— fj;(A?-Al)^ 



UikUklUii , v;^ WifcW'fci + M'ifc'Wfci 



^ 2^2^ (A^ _ A2VA2 _ A21 +2^ 



(2.92b) 



^(A2-A2)(A2-A2) ' ^^ (A2-A2^ 



where the A^'s are the components the p x p diagonal matrix A. The above formulae may 
be checked by verifying that the first evolution equation belonging to the bihamiltonian 
ladder containing 7i2,i can be written as 

L = {L, HoA^^^ = {L , n2,^V'^ = [eu , L] , {n^,M, w) = j ^ V (2.93) 

consistently with (2.75), (2.86) and (2.87). The equation generated by 7^2, « via the second 
Poisson bracket can now be obtained by substituting the expressions for the functional 
derivatives ^^'^ and ^^' into (2.91a,b). The result is somewhat lengthy, so we do not 
display it here. 

We next give the first two Hamiltonians in the bihamiltonian ladder of eq. (2.83); i.e., 

^fc = E(/^^)'^M = r /"trresL^/2= /\^ for A; = l,3, (2.94) 

,^1 /« Jo Jo 
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where the entries of p = diag (pi, . . . , Pp) are ±l's, (see (2.79)). To write down the result, 
it is useful to define the diagonal matrix a by 

a = diag (ai , . . . , CTp) := pA , (2.95) 

and to associate to any p x p matrix v the matrix v given by the formula 

(U),, := j-^ . (2.96) 

[ai + aj) 



Using this notation, we have 



hP = tr a-^w - tr a-^u'^ , (2.97) 



and 



h^ =tr w i w — -u'a — 2au — 2u^ j + tr u i au a + -cr'^u 



(2.98) 



+ tr w^ ( (U^) + 2au + -u a^ 
The first case of equation (2.83) is given explicitly as follows: 

u = {u, 7^^}^2) = {u, n^Y^'^ =au' -2a^u' +[u,u] + [a,w], (2.99a) 

w = {w, 7^^}(2) = {w , n^}^^'^ = aw' - a^u" - uu +\u,w]. (2.99b) 

These equations simplify to the free chiral wave equation in the scalar case p = 1 after 
putting -u = 0, as it should. The analogue of the KdV equation is generated by Tig through 
the second Poisson bracket (2.91a,b). It is straightforward to derive it from the above, but 
the final expression is quite long. 

Remark 2.7. There is no singularity in the formulae (2.92b), (2.96-98) above since the 
diagonal matrix A^ = (— D)"^ has distinct, non-zero eigenvalues. This goes back to 
choosing A = A^. (S> -D as a regular element of the Heisenberg subalgebra at the beginning 
of the construction. In general, for any r, if we continuously deform D to the unit matrix, 
which represents a singular case (type II in the terminology of [7]), then the Hamiltonians 
'Hmr,i in (2.74) (m = 1,2,...) may be expected to become oo and thus disappear from the 
hierarchy. In addition, it follows from a result in [16] that in the D = I case there exists 
only a single rth root of L, up to a scalar multiple. 
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3. Discussion 

This paper was aimed at furttier exploring the integrable systems that can be associ- 
ated to graded regular elements of loop algebras by the method of Drinfeld and Sokolov. 
We have concentrated on the simplest case, given by the nontwisted loop algebra £{gln), 
and have shown that graded regular elements exist only in those Heisenberg subalgebras 
that correspond to partitions into equal blocks n = pr or equal blocks plus a singlet 
n = pr + 1. We further analyzed the first case by taking the grade 1 regular elements and 
proved that the generalized DS reduction results in the matrix version of the r-KdV hier- 
archy of Gelfand-Dickey. We wish to close by mentioning some other interesting models, 
which are related to this KdV type hierarchy in a way that is familiar in the scalar case 
p — 1. 

First of all, one has the modified KdV type hierarchy which is obtained from the 
general AKS system by using the block diagonal gauge. This means that the Hamiltonians 
of the modified hierarchy are still given by the functions of L listed in (2.74), but L (and 
consequently L) is considered to depend on 

C = d + 6 + A, where ^ := diag [^i, . . . , ^^] , (3.1) 

through the Miura map (2.22). The evolution equations of the modified hierarchy are 
generated by means of the free current algebra Poisson bracket carried by the ^^'s, which 
are the basic fields of the modified hierarchy. 

The nonabelian affine (periodic) Toda field equation may be defined as follows. Take 
the basic Toda field to be a G/^i-valued block diagonal matrix g{x,t), 

g = dmg[gi,...,gr\, (3.2) 

where the giS (z = 1, . . . , r) are G/p-valued functions, periodic in the space variable x. In 
addition to the grade 1 regular element A, choose also a grade —1 regular element A and 
define the Toda equation to be the zero curvature equation 

[£+,£_] =0, (3.3a) 

where 

£+ := d+ + g-^d+g + A , £_ := c>_ + g-^Kg , (3.3b) 

and d± = {d^ ± dt). Here, £+ is obtained from the modified KdV operator C in (3.1) 
by substituting 9+ for d and g~^d^g for 9, and therefore we have the same relationship 
between the conservation laws of the present Toda model and the modified KdV hierarchy 

37 



as is familiar in the abelian case [1], [10], [25]. In particular, by the same arguments as in 
[1], we can construct an infinite number of conserved local currents for the Toda equation 
by transforming £_|_ into the Heisenberg subalgebra. Note also that another infinite set 
of conserved currents can be constructed by utilizing the fact that the Toda equation can 
equivalently be written as 

[£+,£_] =0, (3.4a) 

with 

£+ := d+ + gAg-^ , £_ := c>_ - d.g ■ g'^ + A , (3.4b) 

and transforming /}_ into the Heisenberg subalgebra. 

The nonabelian affine Toda model defined above was proposed originally by Mikhailov 
[26] and further studied, e.g., in [27]. More precisely, these authors took A = Ar ^ Ip and 
A = A~^ eg) Ip, where Ip is the p x p unit matrix. These are not regular elements of 
the Heisenberg subalgebra and therefore the DS construction of local conservation laws 
corresponding to graded generators of the Heisenberg subalgebra would not be applicable 
with this choice without modifications. 

A related nonabelian conformal (open) Toda model can be obtained by omitting the 
A dependent terms from C± {C±) in the above. This model is a member of a family of 
models described in [28] and can also be obtained by a Hamiltonian symmetry reduction 
of the Wess-Zumino-Novikov-Witten (WZNW) model. By using the WZNW picture it is 
clear that the W-algebra given by the second matrix Gelfand-Dickey Poisson bracket can 
be realized as the algebra of Noether currents in the nonabelian conformal Toda model. 
Details can be found in refs. [4], [5], [18]. Note that in these papers the group Sin was 
used instead of G/^, but one can impose the 5"/^ constraints tr{D^ui) = 0, tr6' = and 
detg = 1 without changing the essential features of the KdV, modified KdV or Toda 
systems, or their relationship. 

It seems plausible that the nonabelian affine Toda model (as well as an appropriate 
"conformal affine" variant of it) can be viewed as a Hamiltonian symmetry reduction of the 
affine WZNW model, generalizing the abelian case [29]. This Toda model is also known 
to be a reduction of the multicomponent Toda lattice hierarchy [30] . Similarly, the matrix 
r-KdV hierarchy should be related to the multicomponent KP hierarchy [30], [31]. 

The other case in which a graded regular element exists in the Heisenberg subalgebra, 
corresponding to the partition n = pr + 1, has not been pursued in the present work. 
By taking an arbitrary regular element of minimal positive grade it may be verified that 
the generalized DS reduction proposed in [7] leads to a W-algebra which is again equal 
to one of those studied in [4], [5] in the context of WZNW reductions. In these papers a 
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family of W-algebras was associated to the SI2 subalgebras of gin (sin). The W-algebra 
arising from a regular element of minimal grade corresponds to the s/2 subalgebra under 
which the defining representation of gl^ decomposes into p copies of the r-dimensional 
sl{2) irreducible representation plus a singlet. 

Recall that both the s/2 subalgebras of gin and the Heisenberg subalgebras of i{gln) 
are classified by the partitions of n. It is unclear whether there is a general relationship 
between all W-algebras associated to s/2 embeddings and KdV type hierarchies or not, 
since there is a W-algebra for any partition, but graded regular elements exist only in 
exceptional cases. 

The present work was based entirely on the Hamiltonian AKS approach to integrable 
systems. The Grassmannian approach [32] , [24] has also been generalized to other Heisen- 
berg subalgebras than the principal one in [33], [31], [9]. It is clear that the integrable 
systems associated to graded regular elements in arbitrary loop algebras would deserve 
further study from both viewpoints. The starting point could be the explicit description 
of Heisenberg subalgebras recently worked out in [34] . 
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